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In the present article we carry out a molecular dynamics study of the core-softened system and 
show that the existence of the water-like anomalies in this system depends on the trajectory in 
P — p — T space along which the behavior of the system is studied. For example, diffusion and 
structural anomalies are visible along isotherms as a function of density, but disappears along the 
isochores and isobars as a function of temperature. On the other hand, the diffusion anomaly may 
be seen along adiabats as a function of temperature, density and pressure. It should be noted that 
it may be no signature of a particular anomaly along a particular trajectory, but the anomalous 
region for that particular anomaly can be denned when all possible trajectories in the same space are 
examined (for example, signature of diffusion anomaly is evident through the crossing of different 
isochors. However, there is no signature of diffusion anomaly along a particular isochor). We also 
analyze the applicability of the Rosenfeld entropy scaling relations to this system in the regions 
with the water-like anomalies. It is shown that the validity of the Rosenfeld scaling relation for the 
diffusion coefficient also depends on the trajectory in the P — p — T space along which the kinetic 
coefficients and the excess entropy are calculated. 

PACS numbers: 61.20.Gy, 61.20.Ne, 64.60.Kw 



I. INTRODUCTION 



It is well known that some liquids (for example, wa- 
ter, silica, silicon, carbon, phosphorus, and some biolog- 
ical systems) show an anomalous behavior in the vicin- 
ity of their freezing lines |ll421|. The water phase di- 



agrams have regions where a thermal expansion coeffi- 
cient is negative (density anomaly), sclf-diffusivity in- 
creases upon compression (diffusion anomaly), and the 
structural order of the system decreases with increasing 
pressure (structural anomaly) [y, lD] . 

The first anomaly mentioned above is density anomaly. 
It means that density increases upon heating or that 
the thermal expansion coefficient becomes negative. Us- 
ing the thermodynamic relation (dP/dT) v — ap/Kx, 
where ap is a thermal expansion coefficient and Kt is 
the isothermal compressibility and taking into account 
that Kt is always positive and finite for systems in equi- 
librium not at a critical point, we conclude that density 
anomaly corresponds to minimum of the pressure depen- 
dence on temperature along an isochor. This is the most 
convenient indicator of density anomaly in computer sim- 
ulation. 

If we consider a simple liquid (for, example, Lennard- 
Jones liquid), and trace the diffusion along an isotherm 
we find that the diffusion decreases under densification. 
This observation is intuitively clear - if density increases 
the free volume decreases and the particles have less free- 
dom to move. However, some substances have a region in 
density - temperature plane where diffusion grows under 



densification. This is called anomalous diffusion region 
which reflects the contradiction of this behavior with the 
free volume picture described above. This means that dif- 
fusion anomaly involves more complex mechanisms which 
will be discussed below. 

The last anomaly we discuss here is structural anomaly. 
Initially this anomaly was introduced via order param- 
eters characterizing the local order in liquid. However, 
later on the local order was also related to excess entropy 
of the liquid which is defined as the difference between 
the entropy and the ideal gas entropy at the same (p, T) 
point: S ex = S — Sid- In normal liquid excess entropy 
is monotonically decaying function of density along an 
isotherm while in anomalous liquids it demonstrates in- 
creasing in some region. This allows to define the bound- 
aries of structural anomaly at given temperature as min- 
imum and maximum of excess entropy. 

The regions where these anomalies take place form 
nested domains in the density-temperature [6| (or 
pressure-temperature 7]) planes: the density anomaly re- 
gion is inside the diffusion anomaly domain, and both of 
these anomalous regions are inside a broader structurally 
anomalous region. It is reasonable to relate this kind of 
behavior to the orientational anisotropy of the poten- 
tials, however, a number of studies demonstrate water- 
like anomalies in fluids that interact through spherically 
symmetric potentials 22h52|. 

As it was discussed in many works (see, for example, 
the reviews [22j and [47[) the presence of two length scales 
in the core-softening potential may be the origin of water- 



like anomalies: a larger one, associated with the external 
finite repulsion (effective at lower pressures and temper- 
atures), and a smaller one, related to the particle hard 
core (dominant at higher pressures and temperatures). 
In those thermodynamic regimes where the two length 
scales are both partially effective and thus are competing 
with each other, a system of particles interacting through 
such potentials behaves, in many respects, as a mixture of 
two species of different sizes. This leads to the existence 
of two competing local structures: an expanded structure 
characterized by large open spaces between particles, and 
a collapsed structure in which particles are spaced more 
closely. The evolution of these structures under changing 
the thermodynamic conditions can result in the anoma- 
lous behavior. For example, as it was shown in many 
works (see, for example, |54J). the low temperature ther- 
modynamic anomalies of liquid water arises from the in- 
termittent fluctuation between its high density and low 
density forms, consisting largely of 5-coordinated and 4- 
coordinated water molecules, respectively. 

However, it should be noted that in general the ex- 
istence of two length scales is not enough to mark the 
occurrence of the anomalies. For example, for the mod- 
els studied in Ref. [55| it was shown that the existence 
of two distinct repulsive length scales is not a necessary 
condition for the occurrence of anomalous phase behav- 
ior. 

The problem of anomalous behavior of core-softened 
fluids was widely discussed in literature (see, for exam- 
ple, the recent review |22J). It was shown that for some 
systems the anomalies take place while for others do not 
|22J . In this respect the question of criteria of anomalous 
behavior appearance remains the central one. However, 
another important point is still lacking in the literature - 
the behavior of anomalies along different thermodynamic 
trajectories. Here we call as "trajectory" a set of points 
belonging to some path in (P, p, T) space. For example, 
the set of points belonging to the same isotherm we call 
as "isothermal trajectory" or shortly isotherm. 

In our previous work we showed 56, 57] that anomalies 
can exist along some trajectories while along others the 
liquid behaves as a simple one. Taking into account this 
result, it is interesting to study the behavior of the quan- 
tities demonstrating anomalies along the different physi- 
cally significant trajectories (isotherms, isochors, isobars 
and adiabats). This investigation will allow to get deeper 
understanding of the relations between anomalous behav- 
ior and thermodynamic parameters of the system which 
spread light on the connection between thermodynamic, 
structural and dynamic properties of liquids. 



II. SYSTEM AND METHODS 

The simplest form of core-softened potential is the so 
called Repulsive Step Potential which is defined as fol- 



lowing: 



U(r) 




(1) 



where d is the diameter of the hard core, a is the width 
of the repulsive step, and e is its height. In the low- 
temperature limit T = ksT je << 1 the system reduces 
to a hard-sphere system with hard-sphere diameter er, 
whilst in the limit T > > 1 the system reduces to a hard- 
sphere model with a hard-sphere diameter d. For this rea- 
son, melting at high and low temperatures follows simply 
from the hard-sphere melting curve P — cT/a' 3 , where 
c « 12 and a' is the relevant hard-sphere diameter (u and 
d, respectively) . A changeover from the low-T to high-T 
melting behavior should occur for T = 0(1). The precise 
form of the phase diagram depends on the ratio s = cr/d. 
For large enough values of s one should expect to observe 
in the resulting melting curve a maximum that should 
disappear as S —> 1. The phase behavior in the crossover 
region may be very complex, as shown in 48j, [52j . 

In the present simulations we have used a smoothed 
version of the repulsive step potential (Eq. (TTJ)), which 
has the form: 

/ j\ n i 
U(r)=e(-\ +-e(l-tanh(*d(r-a.))), (2) 

where n — 14, fco = 10. We have considered a s = 1.35. 
Here and below we refer to this potential as to Smooth 
Repulsive Shoulder System (SRSS). 

In the remainder of this paper we use the dimensionless 
quantities: f = v/d, P = Pd 3 /e, V = V/Nd 3 = 1/p. As 
we will only use these reduced variables, we omit the 
tildes. 

In Refs. 48j,[52J], phase diagrams of SRSS models were 
reported for er s = 1.15, 1.35, 1.55, 1.8. 

Fig. [Tj shows the phase diagrams that we obtain from 
the free-energy calculations for a s = 1.35. In fact, the 
phase diagrams for a* = 1.15, 1.35, 1.55, 1.8 were already 
reported in Refs. 48j,|52j. We show these phase diagrams 
here too because they provide the "landscape" in which 
possible "water" anomalies should be positioned. 

Fig. Ufa) shows the phase diagram of the system with 
a s = 1.35 in the p — T plane. There is a clear maximum 
in the melting curve at low densities. The phase diagram 
consists of two isostructural FCC domains corresponding 
to close packing of the small and large spheres separated 
by a sequence of structural phase transitions. This phase 
diagram was discussed in detail in our previous publi- 
cations 48j, |52j . It is important to note that there is a 
region of the phase diagram where we have not found any 
stable crystal phase. The results of Ref. [48] suggest that 
a glass transition occurs in this region with vitrification 
temperature T g = 0.079 at p = 0.53. The apparent glass- 
transition temperature is above the melting point of the 
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FIG. 1: (Color online). Phase diagram of the system of par- 
ticles interacting through the potential ((2} with a a — 1.35 in 
p - T (a) and P - T (b) planes. 



low-density FCC and FCT phases. If, indeed, no other 
crystalline phases are stable in this region, the "glassy" 
phase that we observe would be thermodynamically sta- 
ble. This is rather unusual for one-component liquids. 
In simulations, glassy behavior is usually observed in 
metastable mixtures, where crystal nucleation is kinet- 
ically suppressed. One could argue that, in the glassy 
region, the present system behaves like a "quasi-binary" 
mixture of spheres with diameters d and a s and that the 
freezing-point depression is analogous to that expected in 
a binary system with an eutectic point: there are some 
values of the diameter ratio such that crystalline struc- 
tures are strongly unfavorable and the glass phase could 
then be stable even at very low temperatures. The glassy 
behavior in the reentrant liquid disappears at higher tem- 
peratures. 

In the present study a system of particles interacting 



via the potential with "hard" core, repulsive shoulder 
and attractive well is also investigated. This potential 
represents a generalization of our previous SRSS model 
48|,[52( and we call it Smooth Repulsive Shoulder System 
with Attractive Well (SRSS- AW) potential J53|. 
The general form of the potential is written as 

/cr\ 14 

U(r) = el -J +e(A - Aitanh(fci{r-cri}) + 

+ A 2 tanh(fc 2 {r-CT 2 })). (3) 

We consider only the potentials with u\ — 1.35 (see 
Table 1). 

Before to proceed, let us consider the analytic condi- 
tion for core softening proposed in Ref. 58| |: 



and also 



rU"{r) + U'(r) < 0, r 2 > r > n 



17" (r) > 



(4) 



(5) 



for r < n and r > r 2 , where n and r 2 denote two repul- 
sive length scales. In the systems satisfying these condi- 
tions there are two local structures which compete with 
each other. In this case the system behaves as the mix- 
ture of two types of particles with effective radii r\ and r 2 
(see, for example, Ref. [48|). and one can expect to find 
a reentering melting and other anomalous behaviors. 

It may be easily seen that for the potentials (1-3) (see 
Table 1) the conditions (j4]) and (O are satisfied for o\ > 
1.16, so in the case u\ = 1.35 the anomalies do exist 



5J,|53|. 

However, as it was mentioned above, in general the 
conditions (J4j and ([5]) are not enough to mark the occur- 
rence of the anomalies [55| . 

In Refs. 22J, |23( the extensive study of the softness 



dependence of the anomalies for the continuous shoul- 
dered well potential was presented and different criteria 
for the appearance of the anomalies was analyzed. It was 
shown that for the more steeper soft-core the regions of 
the density and diffusion anomalies become more narrow, 
while the region of the structural anomaly is only weekly 
affected, and the disappearance of the density and diffu- 
sion anomalies for the steeper potentials is due to a more 
structured short-range order. At the same time, it may 
be shown that the conditions (j4]) and (El) are satisfied for 
the potentials considered in Refs. 22|,l23(. For the steeper 
potentials the range between r\ and r 2 is more narrow, 
and the effective diameter of the soft core is larger. How- 
ever, as it was shown earlier 52j, the increasing of the soft 
core diameter leads to the disappearance of the anoma- 
lies in the range of the thermodynamic stability of the 
system. It seems that the search of the adequate cri- 
tcrium relating the appearance of the anomalies with the 
intermolecular potential is very interesting and impor- 
tant problem. 
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FIG. 2: (Color online) Family of the potentials with o\ — 1.35 
and different attractive wells. The curves are numerated in 
accordance with Table 1. 



We consider the system with the step a\ = 1.35. The 
parameters k\ = hi = 10.0 are fixed while parameters 
<72, Ao, Ai and A2 are varied to get the different potential 
shape. Three sets of parameters are considered. They 
are summarized in Table 1. Fig. [2] shows the potentials 
for the step width <j\ = 1.35 respectively. The parame- 
ters are chosen in such a way that the depth of attractive 
well becomes larger (see Table 1 and Fig. [2J . Below we 
denote the systems with different parameters as system 
1, system 2 and so on in accordance with the Table 1. 
The choice of the potential parameters was dictated only 
by the convenience for the analysis of the qualitative in- 
fluence of the attraction on the properties of the system. 

In our previous publications [48J, [52J we discussed the 
phase diagrams of several purely repulsive systems, i.e. 
the systems with zero well depth. The complexity of 
these phase diagrams was shown. The systems with at- 
tractive well were considered in [53j. For the complete- 
ness we present the evolution of the system behavior with 
increasing attraction in Figs. [3] and |U 

The transition lines were determined from the free en- 
ergy calculations 5JJ, |60j . 

We simulate the system in NVT ensemble using 
Monte-Carlo method. The number of particles in the 
liquid or gas state simulation was set to 500 or 1000 and 
for crystal phases it varied between 500 and 1000 de- 
pending on the structure. The system was equilibrated 
for 10 6 MC step and the data were collected during 10 5 
MC steps. 
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FIG. 3: (Color online) Phase diagram of the system 2 (Table 
1) in (a) p — T and (b) P — T coordinates. 



In order to find the transition points we carry out the 
free energy calculations for different phases and construct 
a common tangent to them. For the purely repulsive 
potentials we computed the free energy of the liquid by 
integrating the equation of state along an isotherm [59J, 



F(p)-F id {p) 



■Si 



P P{p')-p'k B T 



dp'. In the case of 



potentials which contain an attractive part the situation 
is more complicated because of the possible gas - liquid 
transition. In order to avoid the difficulties connected to 
this transition we carry out calculation of free energies 
at high temperature above the gas - liquid critical point 
and then calculate the free energies by integrating the 



internal energies along an isochor [5 

rT 2 



6C 



F(T 2 )-F(Ti) _ 
k B T 



f^U(T,N,V)d(±). 

Free energies of different crystal phases were deter- 
mined by the method of coupling to the Einstein crystal 

mill. 

To improve the statistics (and to check for internal 
consistency) the free energy of the solid was computed 
at many dozens of different state points and fitted to 
multinomial function. The fitting function we used is 
a P,qT p V q , where T and V — 1/p are the temperature and 
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FIG. 4: (Color online) Phase diagram of the system 3 (Table 
1) in (a) p — T and (b) P — T coordinates. 



specific volume and powers p and q are related through 
p + q < N. The value N we used for the most of calcu- 
lations is 5. 

The diffusion anomaly is also discussed in the article. 
Since the diffusion coefficient can not be measured in 
Monte Carlo simulations, molecular dynamics is applied. 
The core-softened systems are characterized by a com- 
plex energy landscape. This makes difficult to simulate 
the system at low temperatures. In order to avoid this 
problem, parallel tempering technique is used [59j, |60| . 
The diffusion is measured along a set of isochors between 
the densities p = 0.3 and p = 0.8. The temperatures 
used are confined between T — 0.15 and T = 0.8. 

We simulate a system of 864 particles in a cubic box. 
Each parallel tempering run consists of 16 exchanges be- 
tween 8 different temperatures. Between the exchanges 
the system evolves for 4 • 10 6 steps. The first 3 • 10 6 steps 
are used for equilibration. The time step is dt = 0.0005. 
In order to keep the temperature constant Andersen ther- 
mostat is used during the equilibration. Summing up 
all simulations done and taking into account exchange 
of temperatures in parallel tempering runs, we collect 
more then a hundred measurements along each isochor 
which gives good statistics. The diffusion coefficient 



along isochors is approximated by a 9— th order poly- 
nomial of temperature. Then the data are rearranged 
along isotherms. 



III. RESULTS AND DISCUSSION 
Diffusion anomaly 

The diffusion anomaly of the SRSS was discussed in 
several our previous articles. 



52, 5c 



In the Refs. [48 
we showed that the diffusion anomaly takes place at the 
shoulder width o\ = 1.35 while in the work 65| the 
breakdown of the Rosenfeld scaling for this system was 
demonstrated. Further discussion of the Rosenfeld scal- 
ing was reported in the works [56 



571. The main idea 



of these papers is that the anomalous diffusion behav- 
ior present along low-temperature isotherms while along 
isochors diffusion coefficient is monotonous. As a re- 
sult, Rosenfeld scaling is valid along isochors and high- 
temperature isotherms, however, it breaks down for the 
isotherms with anomalies. It means that appearance or 
not of anomalies at some trajectory can cause physically 
different behavior of the system. As a result, following 
different trajectories, we can or can not observe some ef- 
fects. It is of particular importance since experimental 
works deal mostly with isobars and isotherms while theo- 
retical studies - with isochors and isotherms. Taking this 
into account, one can expect that some of the effects ob- 
served along isobars are not visible along isotherms and 
isochors which makes us confused while comparing ex- 
perimental results with theoretical predictions. Here we 
extend the study of anomalies to four different physically 
meaningful trajectories: isotherms, isochors, isobars and 
adiabats and consider both SRSS and SRSS-AW poten- 
tials. 

Fig. EJa) shows the diffusion coefficient of SRSS sys- 
tem along a set of isotherms for ci = 1.35 (potential 1 
in Table 1). One can see that the diffusion coefficient 
demonstrates anomalous behavior for the temperatures 
below T = 0.25. On the other hand, if we look at the 
Fig. [5jb) where the same diffusion coefficient data are 
arranged along isochors as a function of temperature we 
do not observe anomalies - the diffusion is a monotonous 
function of temperature along isochors. However, some 
of the isochors cross. One can see, that the cross of the 
isochors corresponds to the densities between p = 0.4 and 
p = 0.6. Comparing it to the isotherms we see that this is 
the region of anomalous diffusion. It means that even if 
we do not see nonmonotonous behavior of diffusion along 
isochors we can identify the presence of anomaly from 
crossing of the isochors. However, this method seems 
to be technically more difficult since we need to mea- 
sure many points belonging to different isochors rather 
then one isotherm. In Fig. [5jc) the diffusion coefficient 
is shown as a function of pressure along the isochores. 
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FIG. 5: (Color online). Diffusion coefficient of the SRSS sys- FIG. 6: (Color online). Diffusion coefficient of the SRSS-AW 
tern along (a) isotherms and (b) and (c) - isochors. The insert system (system 3 in Table 1) along (a) isotherms and (b) and 
in (a) shows the low temperature isotherms. (c) - isochors. 



One can see that in the range of the densities between 
p = 0.4 and p = 0.6 the slopes of the curves change sign. 
Isothermal and isochoric behavior of diffusion in anoma- 
lous region for SRSS have already been discussed in our 
previous publications [56J, |57j and we give these plots here 
for the sake of completeness. 

If attraction is added to the potential (systems 2 and 3 
in Table 1), no new qualitative features are found, how- 
ever, the anomalies become more pronounced. One can 
see this, for example, in Figs. [5] ((a), (b) and (c)), where 
the diffusion coefficients of the SRSS- AW potential of the 
system 3 are shown along the isotherms and isochores. 

Fig. Ufa) shows the diffusion coefficient along a set of 
isobars as a function of density. The diffusion coefficient 
is again monotonous. The slope of the curves is always 
negative. However, as it can be seen from the insert, the 
slope approaches infinity at low diffusions at pressures 
P = 2.0 and P = 2.5. This corresponds to densities 
0.45 — 0.50 inside the anomalous region. One can imagine 
that if we lower the temperatures along these isobars we 
can observe change of the slope to positive one, however, 
we do not have data for these temperatures. 

Fig. [7[b) shows the diffusion coefficient along isobars 
as a function of temperature. The situation is analo- 
gous to the case of isochors: the curves are monotonous, 
however, they intersect at low temperatures, correspond- 
ing to anomalous region (see the insert in Fig. EJb)). It 
means that if we have the diffusion coefficient along iso- 
bars we can identify the presence of anomalies by mon- 
itoring the intersections of the curves. However, by the 
reasons discussed above this method is not practically 
convenient. 

Figs. [51(a) and (b)) show the diffusion coefficient along 
isobars as functions of density and temperature if attrac- 
tion is added (system 3 in Table 1). As in the case of 
isotherms and isochores, the attraction makes the anoma- 
lies more pronounced. For example, without attraction 
we could find only hints of the diffusion anomaly as a 
function of density (Fig. [7] (a)), however, this anomaly is 
explicit for the system 3 (see Fig. [8J (a)). 

The last physically meaningful trajectory considered 
in the present work is the adiabat. This trajectory 
is defined as constant entropy curve. The entropy is 
calculated as following. We compute excess free en- 
ergy by integrating the equation of states: ^"V ~ 

TT& = kFrfo P P{P V kBT dp'- The excess ^ r °Py 
can be computed via S ex — ~/vF~r" - ^ ne total en- 
tropy is S — S ex + Sid, where the ideal gas entropy is 



^ = fln(T) - \n(p)+\n ^~ h >y ). The last 
term in this expression is constant and is not accounted 
in our calculations. 

The behavior of entropy itself will be discussed below. 
Here we give the diffusion coefficients along the adiabats 
(Fig.El(a)-(c)). 

One can see from the Fig. |H] (a)-(c) that the anomaly 
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FIG. 7: (Color online). Diffusion coefficient of the SRSS sys- 
tem along isobars as a function of (a) density and (b) tem- 
perature. The inserts show anomalous regions in the corre- 
sponding coordinates. 



takes place along adiabats in all possible coordinates 
(look, for example, S = —4.0 adiabat). It means that in 
case of adiabatic trajectory one can identify the anoma- 
lous region monitoring any of three thermodynamic vari- 
ables (P, p, T). However, this trajectory is rather difficult 
to realize in simulation or experiment. 

As one can expect, the attraction does not produce 
qualitative changes in the behavior of the anomalies, but 
makes them more explicit. We do not show the corre- 
sponding figures. 



Density Anomaly 

As we mentioned above density anomaly corresponds 
to appearance of a minima on isochors of the system. The 
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FIG. 8: (Color online). Diffusion coefficient of the SRSS-AW 
system (system 3 in Table 1) along isobars as a function of 
(a) density and (b) temperature. 



isochors of the system with the potential @ (system 1 
in Table 1) are shown in the Fig. [TUTa). It is evident 
from the figure that some of the isochors do demonstrate 
minima. The location of the minimum in the p — T plane 
is shown in Fig. HUT d). 

If we turn to the shape of the isobars themselves for 
the potential © (system 1 in Table 1), that is the de- 
pendence of temperature on density at fixed pressure, we 
do not find any traces of anomalies there (Fig. [TTj) . Like 
in the case of diffusion, the curves have negative slope 
which approaches zero at low temperatures and densities 
corresponding to anomalous regime. However, the curve 
remains monotonous, and it seems that there is no sign of 
density anomaly along isobars. However, using the well 
known thermodynamic relation: 
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FIG. 9: (Color online). Diffusion coefficient of the SRSS sys- 
tem along adiabats as a function of (a) temperature, (b) den- 
sity, and (c) pressure. The inserts show anomalous regions in 
the corresponding coordinates. 
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FIG. 11: (Color online). A set of isobars of SRSS. 
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FIG. 10: (Color online), (a) A set of isochors of SRSS. The 
stars show the location of minimum. The insert enlarges the 
p = 0.45 isochor. (b) The location of the minima on isochors 
in p — T plane. 



N 



t- , where Kt is 



one can find that Kt (f^) v ( §? ) 

the isothermal compressibility. Taking into account that 
Kt is always positive and finite for systems in equilib- 
rium, and using Fig. ITUTa). one can see that ( ^- J > 

if (§y) v < 0. We can conclude that the anomaly does 
exist along isobars for low temperatures, but we do not 
see it in our simulations. 

As it was mentioned above, adding the attraction to 
the potential makes the anomalies more pronounced. To 
illustrate this, in Fig. [12] we show the isobars for the 
system 3 (see Table 1) where the anomalies are clear seen. 

Figs. [13] ((a) and (b)) show adiabats of the SRSS in 
p — T and P — T coordinates. Interestingly, the curves 
seem monotonic, but the slope of the curves at low and 
high densities (pressures) is very different. At the same 
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FIG. 12: (Color online). A set of isobars of SRSS- AW (system 
3 in Table 1). 



time the middle density (pressure) curves demonstrate 
a continues change from the high slope (low-density 
regime) to low slope (high-density regime). However, 

from the well-known thermodynamic relation I 4- ) = 

P 2 ~c~ (w)v an< ^ ^^- [IHlone can see that in the case of 
Fig. [|2] (a) the anomaly does exists, but it is not seen 
because of the insufficient accuracy of calculations along 
the adiabats. On the other hand, due to the relation 



\dPJ 



— (^K) r, there is also anomalous behavior in 

cp V oT ) P 

Fig. \m (b). 

As one can expect (we do not represent these figures 
here for the sake of brevity), the anomalies are much 
better seen for the systems with the attractive potentials 
(systems 2 and 3 in Table 1). 
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FIG. 13: (Color online). Adiabats of the SRSS in (a) p - T 
and (b) P ~T coordinates. 



FIG. 14: (Color online). Excess entropy of SRSS fluid along 
(a) isotherms and (b) isochors. 



Structural Anomaly 

Structural anomaly region can be bounded by using 
the local order parameters or by excess entropy minimum 
and maximum. Here we apply the definition via excess 
entropy. 

The behavior of excess entropy is qualitatively anal- 
ogous to the behavior of diffusion coefficient. Because 
of this we briefly describe it here noting that most of 
the conclusions about diffusion coefficient along different 
trajectories can be applied to excess entropy as well. 

Figs. Q3] ((a) and (b)) show the excess entropy along 
isotherms and isochors for the purely repulsive poten- 
tial ([2]). Like for the diffusion coefficient, excess entropy 
demonstrates anomalous grows in some density range 
at low temperatures. At the same time excess entropy 
is monotonous along isochors. However, the curves for 
different isochors cross which indicates the presence of 



anomaly. 

The excess entropy along isobars is monotonically de- 
creasing function of density and monotonically increasing 
function of temperature (see Figs. [Pol fa) and (b))). How- 
ever, the curves cross at low temperatures indicating the 
presence of anomalies as it was discussed for the case of 
diffusion. As in the case of diffusion anomaly, the struc- 
tural anomaly is not seen along the isobars for the purely 
repulsive potential (system 1 in Table 1), however, one 
can expect that in our simulation we could not reach the 
anomalous region. If the attraction is added to the po- 
tential, the anomaly becomes explicit even as a function 
of density (see Figs. IWfa) and (b))). 

It is important to note that the range of densities of 
structural anomalies is wider then the one of diffusion and 
density anomalies which is consistent with the literature 
data for core-softened systems [y, |7(. Figs. IrTT(a). (b) 
and (c)) represent the locations of the anomalies lines on 
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FIG. 15: (Color online). Excess entropy of SRSS fluid along 
isobars as a function of (a) density and (b) temperature. The 
insert in (b) shows the cross of the curves at low temperatures. 



FIG. 16: (Color online). Excess entropy of SRSS- AW fluid 
along isobars as a function of (a) density and (b) temperature 
(system 3 in Table 1). 



the phase diagrams of the systems in Table 1 in p — T 
plane. 



IV. ROSENFELD SCALING 

In 1977 Rosenfeld proposed a connection between ther- 
modynamic and dynamical properties of liquids [6l|, [62 1 . 
The main Rosenfeld's statement claims that the trans- 
port coefficients are exponential functions of the ex- 
cess entropy. In order to write the exponential rela- 
tions Rosenfeld introduced reduction of the transport 
coefficients by some macroscopic parameters of the sys- 
tem 
D* 



For the case of diffusion coefficient one writes: 

1/3 

D — £,.__m/, , where m is the mass of the parti- 



(t B T/m)'" 

cles. The Rosenfeld scaling rule can be written as: 



D* = A 



^BS C 



(7) 



where A and B are constants. 

In his original works Rosenfeld considered hard 
spheres, soft spheres, Lennard- Jones system and one- 



component plasma [61|, |62j. After that the excess en- 
tropy scaling was applied to many different systems in- 
cluding core-softened liqui ds Ila.[^l63l - l65| . liquid metals 
66J |67| , binary mixture s 1681 l69il , ionic liq uids 7Cl [71 1 , 
network- forming liquids 63J, [7fJ , water [72| , chain fluids 



we showed that 



[73j | and bounded potentials 
In our recent publication 
for the case of the core-softened fluids the applicability of 
Rosenfeld relation depends on the trajectory. In partic- 
ular, Rosenfeld relation is applicable along isochors, but 
it is not applicable along isotherms. 
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The breakdown of the Rosenfeld relation along 
isotherms can be seen from the following speculation. 
The regions of different anomalies do not coincide with 
each other. In particular, in the case of core-softened 
fluids the diffusion anomaly region is located inside the 
structural anomaly one. It means that there are some re- 
gions where the diffusion is still normal while the excess 
entropy is already anomalous. But this kind of behavior 
can not be consistent with the Rosenfeld scaling law. 

However, from this speculation it follows that the 
Rosenfeld scaling should hold true along the trajectories 
which do not contain anomalies, i.e. isochors and isobars. 



In our recent publication 57| we considered the Rosen- 
feld relation along isotherms and isochors. Here we bring 
these trajectories for the sake of completeness and add 
the verification of the Rosenfeld relation along isobars 
for the purely repulsive potential ^ (Figs. UST a) - (c)). 
One can see that the Rosenfeld relation does break down 
along isotherms which is consistent with the speculation 
above. At the same time it holds true along both isochors 
and isobars which is consistent with the monotonous be- 
havior of both diffusion coefficient and excess entropy 
along these trajectories. 

In Figs. ITW /a) - (c)) we show the Rosenfeld relation 
for the system 3. One can see that the behavior is the 
same as in the case of the purely repulsive potential. 

This observation makes evident that Rosenfeld relation 
is valid only along the trajectories without anomalous 
behavior. 



V. SUMMARY AND DISCUSSION 

To summarize, in the present article we carry out a 
molecular dynamics study of the core-softened systems 
(SRSS and SRSS-AW) and show that the anomalous be- 
havior can be seen only along some particular trajectories 
in (P, p, T) space along which the behavior of the system 
is studied. For example, diffusion and structural anoma- 
lies are visible along isotherms as a function of density, 
but disappears along the isochores and isobars as a func- 
tion of temperature. On the other hand, the diffusion 
anomaly may be seen along adiabats as a function of 
temperature, density and pressure. Density anomaly ex- 
ists along isochors, isobars and adiabats. However, if a 
single curve does not demonstrate the anomalous behav- 
ior, having a set of the curves, one can see the presence 
of anomalies via the curves crossing. 

We also analyze the applicability of the Rosenfeld en- 
tropy scaling relations to this system in the regions with 
the water-like anomalies. It is shown that the validity of 
the Rosenfeld scaling relation for the diffusion coefficient 
also depends on the trajectory in the P—p—T space along 
which the kinetic coefficients and the excess entropy are 
calculated. In particular, it is valid along isochors and 
isobars, but it breaks down along isotherms. The break- 
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FIG. 18: (Color online). Rosenfeld relation for SRSS along 
(a) isochors, (b) isotherms, and (c) isobars. 



FIG. 19: (Color online). Rosenfeld relation for SRSS- AW 
along (a) isochors, (b) isotherms, and (c) isobars. 
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down of the Rosenfeld relation along isotherms is related 
to the fact that the boundaries of different anomalies do 
not coincide with each other. 

The influence of attraction on the diffusion anomaly 
is discussed. It is shown that the attraction makes the 
anomalies much more clear and may be found even in the 
cases when the anomalies are hardly seen for the purely 
repulsive potential due to simulation limitations. 

The question is whether this type of behavior is uni- 
versal for all types of systems demonstrating the anoma- 
lous behavior. It is now widely believed that the water 
anomalies are related with the hypothesized liquid-liquid 
critical point, the terminal point of a line of first-order 
liquid- liquid phase transitions [8|, |9|, l76l478| . The anoma- 
lies arise from crossing the Widom line emanating from 
the hypothesized liquid-liquid critical point (LLCP) [771 - 
l79j |. In particular, it was shown [77J that the dynami- 
cal crossover from Arrhenius to non-Arrhenius behavior 
(strong-fragile transition) takes place both for water and 
Jagla [33] model. The Jagla model with attraction dis- 
plays (without the need to supercool) a liquid-liquid co- 
existence line that, unlike water, has a positive slope [77]. 
We believe that our model should show the similar phase 
behavior as the Jagla one, however, with the LLCP in 
the deeply supercooled region [53j. 

For example, in Figs. [30] (( a ) - (c)) we show the be- 
havior of the isothermal compressibility Kt, the isobaric 
heat capacity Cp, and the thermal expansion coefficient 
ap for the system 1 in Table 1. One can see that the 
system demonstrates the behavior compatible with hy- 
pothesis of the liquid-liquid critical point in the deep su- 
percooled region of the phase diagram below the homo- 
geneous nucleation line (compare, for example, with Fig. 
1 in Ref. [79( ) . In Fig. [21] one can see the dynami- 
cal crossover from Arrhenius to non-Arrhenius behavior. 
At high temperature, D exhibits an Arrhenius behavior, 
whereas at low temperature it follows a non-Arrhenius 
one. From Figs. [20] and [21] one can easily see that the 
anomalous lines correspond to the regions of anomalies 
in Fig. [T] (b). This result is consistent with the case of 
the Jagla model [z3]( see Fig. 5 (c) in Ref. (z3|), so we 
believe that in our case the Widom line should have the 
same slope, as in this model, and we can expect that the 
anomalous behavior discussed in this article will be the 
same for the Jagla model too, however, the case of water 
needs an additional investigation. 
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FIG. 20: (Color online). Isothermal compressibility Kt (a), 
isobaric heat capacity Cp (b), and thermal expansion coef- 
ficient ap (c) as a function of temperature along isobars for 
the system 1 in Table 1. 
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FIG. 21: (Color online). Diffusivity D along isochores for the 
system 1 in Table 1. 
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